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Abstract 

A leading order determination of the gluon density in the proton has been per- 
formed in the fractional momentum range 1.9 • 10~ 3 < x g/p < 0.18 by measuring 
multi-jet events from boson-gluon fusion in deep-inelastic scattering with the HI 
detector at the electron-proton collider HERA. This direct determination of the 
gluon density was performed in a kinematic region previously not accessible. The 
data show a considerable increase of the gluon density with decreasing fractional 
momenta of the gluons. 
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1. Introduction 



Deep-inelastic lepton-nucleon scattering experiments have played a fundamental role in 
reaching an understanding of the structure of matter. Ever since the discovery of the 
proton's parton content in the late 1960s, extensive studies have been made at accelerators 
providing increasingly higher energies to obtain more detailed knowledge of the parton 
properties within nucleons. Although it was soon realized that about 50% of the nucleon 
momentum was carried by gluons, a direct measurement of their momentum distribution 
has so far been restricted to large fractional momenta. Instead information about the 
gluons has been extracted from measurements of the sea-quark distribution. This relies 
on the assumption that sea-quark pairs are produced by quark-antiquark production from 
gluons in an evolution process. 




Figure 1: Generic Feynman diagrams for a) the boson-gluon fusion process and, b) the 
QCD-Compton process 

The electron-proton collider HERA has considerably extended the kinematic region 
available to investigations of nucleon constituents. In particular it permits measurements 
of processes directly initiated by gluons in the proton. In the boson-gluon fusion process 
(Fig. la) a gluon from the proton interacts with the virtual boson from the electron to 
produce a quark-antiquark pair. In the range of momentum transfer (Q 2 ) considered 
in this analysis photon exchange is dominant and the contribution from Z° exchange 
can be neglected. Typically the final state contains two jets in addition to the proton 
fragment which, to a large extent, disappears undetected down the beam pipe. Such 
events are denoted as (2+1) jet events. The hard sub-system is the system produced by 
the interaction of the virtual photon with a parton in the proton. The (2+1) jet cross 
section also receives a contribution from the QCD-Compton process (Fig. lb) where a 
gluon is emitted from the interacting quark. The cross section of (2+1) jet events can in 
leading order schematically be written as 

o- 2+ i oc a s (A ■ g + B ■ q) 

where a s is the strong coupling constant and g and q stand for the gluon and quark 
densities. The coefficients A and B can be calculated in pertubative QCD. The (2+1) 
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jet rate at HERA offers a possibility to determine the variation of the strong coupling 
constant with the momentum transfer of the exchanged boson. Such a study has been 
carried out ]l[ based on the assumption that the quark and gluon densities in the proton 
are sufficiently constrained at lower momentum transfer by previous experiments. The 
analysis of the gluon density presented here, on the other hand, relies upon a reasonably 
good knowledge about the strong coupling constant and the quark densities. In the 
kinematic region used for this analysis, the multi-jet processes are initiated by partons 
carrying a momentum fraction of the proton of the order of 10~ 2 , a region in which the 
quark densities have been accurately measured in previous experiments [||, [§. 

2. The HI Detector 

A detailed description of the HI detector can be found elsewhere J4[. Closest to the 
interaction point there are tracking devices surrounded by a calorimeter consisting of an 
electromagnetic and a hadronic section. Outside these detectors a superconducting coil 
provides a magnetic field parallel to the beam line and, finally, the instrumented magnet 
iron gives a rough measurement of the energy leaking out of the calorimeter, and signals 
the presence of a muon track. Here we only give a brief description of the detector parts 
which are of relevance to the measurement of the gluon density. 

The backward electromagnetic calorimeter (BEMC) covers the angular range 155° < 
9 < 176° where 9 is defined with respect to the proton beam direction, which is in the 
following called the forward direction. Making use of the track information from the 
backward proportional wire chamber (BPC), located right in front of the BEMC, and the 
reconstructed vertex position, the energy and scattering angle of the scattered electron 
can be measured in the kinematic region considered here. The BEMC consists of a 22.5 
radiation length deep lead-scintillator sandwich stack, each layer read out by two pairs 
of oppositely positioned wavelength shifter bars. This system gives an energy resolution 
of a(E)/E « 0.1/y/E[GeV] © 0.42/£[GeV] © 0.03. By adjusting the measured electron 
energy spectrum to the kinematic peak, the BEMC energy calibration is known to an 
accuracy of 1.7% |5[]. The BPC has four wire planes, giving a spatial resolution of 1.5 
mm, which together with the precision in the vertex reconstruction, results in a resolution 
of the polar angle of 2.5 mrad. 

For the measurement of the hadronic final state the liquid argon calorimeter was 
used, which extends over the angular range 4° < 9 < 153° with a complete azimuthal 
coverage. The total depth varies between 4.5 and 8 interaction lengths. Measurements 
on test beams have given a hadronic energy resolution of a(E)/E « 0.5/ ^ E[GeV] ©0.02 
|J. The absolute hadronic energy scale has been measured to a precision of 5% from 
the transverse momentum balance between hadronic jets and the scattered electron in 
deep-inelastic scattering (DIS) events. 
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3. Kinematics 



In 1993 HERA was operated to collide 26.7 GeV electrons with 820 GeV protons, resulting 
in a center-of-mass energy of 296 GeV. Fig. 1 shows Feynman diagrams of the boson- 
gluon fusion process (BGF) and the QCD-Compton process (QCD-C) with the relevant 
kinematic variables indicated. 

The overall kinematics of an event can be determined from two independent Lorentz 
invariant variables. This can be any two of the Bj0rken scaling variables x and y, the 
photon momentum transfer squared Q 2 , and the invariant mass squared of the hadronic 
system W 2 . These variables are defined in terms of the four- momenta of the incoming 
proton, P, the incoming and outgoing electron, p e and p' e , and the exchanged photon, q. 
Experimentally they are deduced from measurements of the energy, E' e , and polar angle, 
9 e , of the scattered electron according to the following relations (neglecting the electron 
and proton masses): 

Q 2 = -q 2 = -Oe-p'J 2 = 4E e E' e cos 2 (9 e /2) 
y=^ = l-{E> e /E e )sm\e e /2) 

P-Pe 



2P ■ q ys \ x 

The center-of-mass energy squared is given by s = 4E e E p , where E e and E p are the 
energies of the incoming electron and proton, respectively. Note that Bj0rken-x is not 
identical to the fractional momentum Xi/ P of the parton i initiating the hard sub-process, 
with i standing for either a quark (q) or a gluon (g). Xi/ P is related to the invariant mass 
squared of the hard sub-system (s) in the following way: 

s + Q 2 / § 



X i/P = — = x I 1 + Tyi J C 1 ) 



ys V. Q 

where s can be measured either directly from the invariant mass of all particles Pj of the 
two jets belonging to the hard sub-system 



52 pj 



(2) 



or from the jet directions in the hadronic center-of-mass system, which is the rest system 
of the exchanged photon and the proton: 

§ = W^e-Wm) (3) 

where rjl and 77! are the parton pseudo-rapidities of the hard sub-process[] and rj = 
— lntan(6'/2). Experimentally the parton directions and momenta are approximated by 
the directions and energies of the jets produced in the fragmentation process. 

1 Variables v in the hadronic center of mass system are denoted as v* while variables in the photon- 
parton center of mass frame are denoted as v. The direction of the z-axis in the hadronic center of mass 
system is defined as the direction of the exchanged photon. Note that the transverse momentum = fir. 
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4. Trigger and Data Selection 



The analysis is performed with deep-inelastic scattering (DIS) events which are selected 
in the detector by applying the trigger requirement for an electron of a local energy 
deposition (cluster) of more than 4 GeV in the BEMC detector. The total event sample 
used in this analysis corresponds to an integrated luminosity of 242 nb _1 from the data 
sample collected in 1993. After the event reconstruction has been performed, a clean DIS 
sample is obtained by applying the following requirements: 

• The scattered electron energy had to satisfy E' e > lOGeV which corresponds to 
y < 0.625. This removes most of the background due to the photoproduction reactions 
where an electromagnetic shower in the hadronic final state is misidentified as an electron 
from the DIS process. 

• The track coordinates in the BPC had to match the cluster center-of-gravity in the 
BEMC to within 4 cm and the lateral spread of the cluster had to be less than 4 cm in 
radius. 

• The momentum transfer squared had to be in the range 12.5 < Q 2 < 80 GeV 2 which 
ensures that the scattered electron is well contained in the BEMC. 

• In order to get an adequate resolution in y from the measurement of the scattered 
electron it was required that y > 0.05. 

• At least one charged track from the hadronic final state was demanded for the 
determination of the vertex position along the beam, which had to be within 30 cm from 
the nominal position in order to suppress beam-induced background. 

After these cuts and the jet selection described below no indication was observed of 
remaining background from photoproduction in our data. This was confirmed by a Monte 
Carlo simulation of photoproduction processes. 

Efficiencies for passing the above cuts have been determined from data using the 
redundancy in our apparatus. Further details can be found in ||. The 4% error in the 
combined efficiency will contribute to the overall normalization error. 

The (2+1) jet selection is carried out in two steps. In the first step the jet find- 
ing is performed by applying a cone algorithm H on calorimeter clusters with AR = 
a/ Ai] 2 + A$ 2 = 1 first boosting to the hadronic center-of-mass system. A77 and A$ are 
the differences in pseudo-rapidity and azimuthal angle between two cluster pairs respec- 
tively. Exactly two jets with a transverse jet energy p^ > 3.5GeV are required. The 
hadronic center-of-mass system is a natural frame for a QCD analysis as it is not sensitive 
to the pt of the hadronic system due to the transverse momentum balance of the electron 
in the laboratory system. It is also a preferred frame for the px based cone algorithm. 
The ratio of (3+1) to (2+1) jet events is less than 15%. 

In the second step following cuts are imposed on the jet system of the event: 

• Both reconstructed jets have to fall inside the angular range 10° < 8j et < 150° in the 
laboratory system so as to be within the volume covered by the liquid argon calorimeter. 
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A further purpose of the lower cut in angle is to remove the very forward region which is 
dominated by the proton fragments and initial state parton radiation which may otherwise 
give rise to a separate jet. 

• The difference in pseudo-rapidity between the two jets in the laboratory frame must 
satisfy A77 < 2. This corresponds to a cut on the scattering angle in the photon-parton 
center of mass system of \fj\ < 1 or 40° < 9 < 140° in the kinematic range used for this 
analysis. Note that this cut, except for the highest x g / p covered by this analysis, is more 
restrictive than the angular cut discussed above. 

• The corrected invariant mass (see section 6) of the two jets of the hard scattering 
system is required to be y/1 > 10 GeV to ensure well defined jet structures (see Fig. 2). 
This cut, together with the jet angular cut Arj < 2, implies p^> > 3.2 GeV and is therefore 
well matched to the p^ cut of 3.5 GeV in the application of the cone algorithm. In addition 
both methods are used to reconstruct s, as given by equations (2) and (3). It is required 
that the difference in |Avi| for the two methods is smaller than 10 GeV. 

With these cuts a sample of 328 (2+1) jet events covering the fractional momentum 
range 0.002 < x g / p < 0.2 and the Bj0rken-x range 0.0003 < x < .0015 is obtained. In the 
kinematic region which was chosen for this analysis, according to the cuts defined above, 
the ratio of BGF and QCD-Compton cross sections is approximately 3/1, independent of 

5. Monte Carlo Generation 

In order to calculate cross sections from matrix elements and to correct for parton shower, 
hadronisation, and detector effects, the Monte Carlo event generator LEPTO 6.1 || was 
used, with the parton density function MRS H [f|. Although it would have been more 
natural for this analysis to use leading order parametrisations of the density functions it 
will be explained in the discussion of errors (section 7) that this is not essential. The 
strong coupling constant a s was calculated in first order QCD with Aqcd = 250 MeV for 
5 flavours. It was checked that LEPTO delivers cross sections which are consistent with 
leading order cross section predictions by the PROJET JTD[] and DISJET [[LI]] programs 
within 5% in the kinematic region covered by this analysis. 

The LEPTO generator is based on QCD matrix element calculations up to first order in 
a s , and the inclusion of parton showers accounts for higher order processes (MEPS). The 
QCD matrix element gives divergences for soft and collinear emission which is technically 
avoided by defining a smallest invariant mass rriij between any two partons, including the 
remnant. The parametrisation used for the cut in is always 2 GeV above the region 
in phase space where the 2+1 jet cross section would exceed the total cross section. With 
such a cut, matrix element calculations are used to generate events over the maximum 
possible phase space, thereby giving access to regions of small x g / p . The Monte Carlo 
matrix element cut-off is always significantly below 10 GeV which is the constant mass 
cut applied in the event selection. Therefore the full phase space used in this analysis is 



8 



-2 



2 



-2 2 4 



AO 



At) 



t3 q — 



g 
> 

o 



c) 



r+t 





AO 




At) 



5 9 

a 

£ 8 

O 7 

§ 6 

> 

5 

1 4 



e) 





AO 

Figure 2: TTie transverse energy flow in the laboratory frame as a function of a) A$ and 
b) At] for the most backward going jet, and as a function of c) A$ and d) A77 for the 
most forward going jet, using the respective jet axis as reference. The energy flow in the 
region 2 < 77 < 3 is shown with the jet axis of the most forward going jet as reference 
in e). The sense of $ is away from the scattered electron for the jet at the smallest 77. 
The points represent the (2+1) jet data sample and the errors are statistical whereas the 
histograms give the predictions of the MEPS model. 
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covered by the matrix element. There is, however, some migration of events with invariant 
masses below 10 GeV. The contribution from the zeroth order process (QPM events) is 
about 60% at the generator level. In the final event sample a 2% contribution from these 
processes is observed, generated as (1+1) jet events, but wrongly reconstructed as (2+1) 
jet events. This indicates that the contribution from the phase space not covered by the 
matrix element is small. 

Although the matrix element calculation used in LEPTO assumes massless quarks, 
the program contains a parametrisation which takes the heavy quark mass thresholds 
into account. In the kinematic region of this analysis, this effect reduces the LEPTO 
BGF cross section by 11% on the average. 

A generated sample corresponding to about four times the statistics of the data was 
subject to full detector simulation, followed by event reconstruction. 

A number of control plots have been made to compare the predictions of the Monte 
Carlo program with the experimental data. Fig. |2]a-d show the transverse energy flow 
as a function of the azimuthal angle, and pseudo-rapidity, rj, for the two jets using 
the respective jet axis as reference. Particular attention was paid to the Monte Carlo 
description of the transverse energy flow in the forward region, 2 < rj < 3, which might 
contain a large contribution from initial state QCD radiation. Although we have previ- 



ously observed (E| that the MEPS model is not able to reproduce the energy flow in this 
rapidity region for an inclusive DIS sample, it is demonstrated in Fig. |2]e that the energy 
flow of the selected (2+1) jet sample is well described by the MEPS model. It can be 
noted that not only the jet profile exhibits good agreement between data and the Monte 
Carlo sample but also the level of the underlying energy flow. Other control plots, such as 
the rapidity distributions of the two hard jets and their transverse energy spectra, as well 
as the Q 2 , W 2 , y and x distributions, give further evidence for the ability of the MEPS 
model to reproduce the data. 

In order to estimate the model dependence of the BGF acceptance LEPTO MEPS was 
compared with other models available. It was found that the MEPS-, CDM 4.03- (colour 
dipole model) and HERWIG 5.7 [14[] models all give acceptances for BGF-events 



agreeing within 15% at the hadron level. The CDM model also gives a good description 
of the jet profiles but it is not able to give a satisfactory description of the cross section 
as a function of kinematic variables and jet variables as previously observed with less 



statistics in ||15|| . The discrepancies originate from non-BGF processes. In HERWIG 
there is no matching between the zeroth and first order processes, preventing an estimate 
of the migration of background. 



6. Reconstruction of x g / p and Unfolding of the Gluon Density 

The Monte Carlo sample was used to investigate the correlation between x g / p of the gluon 
calculated from the hard partons originating from the QCD matrix element and x r g e ,^ 
calculated from the jets measured in the detector. It is essential to use a jet reconstruction 
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Figure 3: The relative error in the reconstruction of the fractional gluon momentum, 



Axg/p/x g / p , where Ax g / P 



x 



g/p 



— Xg/p, for BGF events generated by the MEPS model. 



In a) equation has been used, in b) equation (|5[) and in c) the combination of the two. 
The correlation between the reconstructed and true x g / p for the combined method is shown 
in d). 



algorithm which provides good separation between the spectator jet and the jets of the 
hard sub-system. The event sample populates the region of small x and Q 2 values, the 
spectator jet will carry away most of the available energy. A misassignment of one of 
the relatively energetic particles from the spectator jet into the hard sub-system has a 
strong impact on the reconstructed s value and will distort the x g / p reconstruction. It is 
therefore important to optimize the resolution parameters of the jet algorithms to give 
the best possible performance in this respect. 

The jet definition and detector effects lead to systematic shifts of the order of 30% 
in the determination of x g / p and s (see equation [I]). The shifts have been determined 
separately for the two methods of extracting x g / p (equations (H) and (^)). Figs ||a and b 
show the relative error in the x g / p reconstruction after correcting for the shift. 
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Comparing the results of the two methods it can be seen that the fitted Gaussian 
distributions result in a superior resolution for the method based on equation (||) which, 
however, exhibits a more pronounced tail. 

For properly reconstructed events both methods are expected to give consistent results, 
whereas a misassignment of particles by the jet algorithm might have different impacts 
on the s (and therefore x g / p ) reconstruction from the two methods. Also for (1+1) jet 
events, where an additional jet arises from high px initial parton emission, the methods 
do not give the same values of s. Therefore an improved result is obtained by selecting 
only events for which the result of the two methods agree within the resolution. It has 
been demanded that the absolute value of the difference between the y/l values extracted 
from equations (0) and (0) is |A\/I| < 10 GeV and then simply taken the mean value of 
the two reconstructed values to give the combined x g / p . Due to this cut about 20% of the 
events are removed. 

As shown in Fig. ^p, the result of the combined method shows a considerably improved 
resolution compared with the methods based on equations (g) and (0) separately. In 
Fig. |]d the correlation between the reconstructed and true x g / p is given for the combined 
method. From here on the combined method is used to extract x g / p . 

The observed (2+1) jet cross section as a function of x T g ? p is interpreted as the sum of 
a gluon initiated (BGF) and a quark initiated (QCD-C plus QPM) part, according to 

a 1bsX x ljp) = j * M(x g / p ,x r g Y p )xg/ p g(x g / p )dx g /p + C (x r q f p ) + \x r q f p ) 

The quark initiated contributions represent the background obtained from the LEPTO 
Monte Carlo program. The integral represents the BGF part and expresses the convolu- 
tion of the gluon density x g / p g(x g / p ) with the function M which besides the measurement 
process describes effects due to the QCD matrix element, parton showering and hadroni- 
sation. The unfolding procedure with regularisation described in flTj| is used to determine 



the gluon density. The ge covered by the kinematic region is subdivided into five 

bins. The bin boundaries (see Tab. 1) were optimized in order to minimize the correlations 
between the bins. 











diffr active 


x 5 /p-bin range 




< Q 2 > [GeV 2 ] 


<Pt> [GeV] 


contribution [%] 


0.0019 - 0.0061 


9.04 + 2.95 


23 


5.0 


16 + 9 


0.0061 - 0.012 


7.40 ± 1.28 


29 


5.4 


10 + 5 


0.012 - 0.030 


4.03 + 1.06 


30 


6.0 


6 + 3 


0.030 - 0.052 


1.42 + 0.60 


35 


6.6 


9 + 4 


0.052 - 0.18 


1.18 + 0.80 


36 


9.8 


9 + 6 



Table 1: Results and kinematic characteristics. Errors include statistical and systematic 
uncertainties, excluding a 11% global normalization uncertainty. 
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Figure 4: The measured gluon density at an average Q 2 of 30 GeV 2 as a function of 
the fractional gluon momentum. The error bars reflect the statistical errors and the total 
errors respectively. Not included is the global normalization uncertainty of 11%. The solid 
curve shows a fit to the data points as explained in the text. 



The results for the gluon density as a function of the fractional gluon momentum are 
presented in Fig. 4. The inner and outer error bars on the data points reflect the diagonal 
elements of the statistical and full covariance matrix respectively. Not shown is the 11% 
global normalization uncertainty. Additional information is given in Table 1. Results 
are quoted here at a scale of 30 GeV 2 corresponding to the average Q 2 value of the data 
sample. Another possible choice for the QCD-scale would be the average p*^ of the jets 
in the hadronic center-of-mass system, which for the present data sample is 40 GeV 2 . 
Comparing the GRV gluon parametrisation |L7j (see below) for Q 2 values of 20 and 40 
GeV 2 gives less than 10% variation in our range of x g / p . 

A considerable rise in the gluon content of the proton x g / p g(x g / p ) with decreasing 
fractional momentum of the gluon is observed. In Fig. 4 a fit of the data points to 
the parametrisation x g / p g{x g / p ) = A g x g f(l — x g / p ) 9 is shown taking the full covariance 
matrix into account (see Tab. 2). The parameter C g was kept fixed at the value of 5. 
The values of the other parameters as obtained from the fit are A g = 0.33 ± 0.16 and 
B g = —0.63 ± 0.12 with x 2 / n f — 0.99. A value of C g = 9 lead to variations of A g and B g 
within the stated errors. 
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7. Discussion of Errors 



The unfolding yields a full covariance matrix of statistical errors which includes informa- 
tion on correlations between bins due to migration effects. In this analysis, the correlations 
are small (see Tab. 2), indicating adequate resolution for choosing five bins in x g / p . For 
a quantitative estimate of systematic uncertainties we varied several parameters and re- 
peated the full analysis including the unfolding step. For each variation, a covariance 
matrix was calculated. The total systematic covariance matrix is the sum of the matrices 
of the individual variations. The following effects were studied. 

The absolute hadronic energy calibration (affecting the measurement of ]f T and s) 
was varied by its present precision of 5%. We also varied the electron energy calibration 
(affecting the measurement of the Bj0rken variable y and the total hadronic invariant 
mass W) by its measured precision of ±1.7%. 

The jet transverse momentum cut in the cone algorithm used for jet finding was varied 
between 3 and 4 GeV. The jet angular cut A77 nominally 2 was varied between 1.5 and 
2.5. 

In section 5 it was shown that the MEPS model gives an excellent description of the 
data sample used in this analysis. However, comparing the distribution of the absolute 
difference of ^/I reconstructed according to equations (2) and (3) (|Av£|, see section 6) 
an excess of data events is observed compared to the MEPS model at values of |Aa/^| 
larger than 15 GeV. Below this value Monte Carlo and data show good agreement. The 
excess in data is due to events with at least one jet in the forward region, a region where 
the Monte Carlo description of data has been shown to be problematic |L2|. However, 



even for the events rejected by the cut |A\/I| < 10 GeV the MEPS model gives a good 
description of the energy flows. For systematic studies, the |A\/I| cut is varied between 
7 and 15 GeV. 

QCD-Compton cross sections using two leading order (LO) quark density parametri- 
sations (GRV [|17j and CTEQ3L [0) have been compared with the next to leading order 
(NLO) MRS H function used in the MC model. The discrepancies are largest at the low 
x q /p but do not exceed 15% in the kinematic domain used. As a conservative estimate of 
the systematic uncertainty related to the subtraction of the QCD-Compton contribution 
a variation of ±25% was used to determine the contribution to the systematic error. 



( 6.220 

0.087 0.978 

-0.257 -0.064 

0.033 -0.065 

\ 0.032 0.013 



0.293 

-0.017 0.214 
-0.049 0.003 



\ 



0.366 / 



( 8.670 

0.858 1.632 

0.976 0.584 

0.243 0.038 

\ -0.407 0.108 



1.117 

0.050 0.354 
-0.162 0.044 



\ 



0.642 / 



Table 2: Covariance matrix of statistical errors (left) and total (sum of statistical and 
systematic) errors (right). Not included is a global normalization uncertainty of 11%. 
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Tab. 2 shows the total covariance matrix which is the sum of the statistical and 
systematic covariances. The relative systematic errors vary between 17% in the lowest 
bin of Xg/p and 44% in the highest bin. The dominating contribution to the systematic 
errors in the lowest x q / p bin is 13% coming from the uncertainty in the hadronic energy 
calibration of the liquid argon calorimeter, while the dominating systematic error in the 
highest x q /p bin is 35% due to the variation of the |A\/I| cut. 

In addition there is a global normalization uncertainty of 11% which arises from the 
following sources: a 4.5% uncertainty in the luminosity measurement |19[ directly propa- 
gates into the gluon density measurement, as does a 4% uncertainty in various detector 
efficiencies mentioned in section 4. Finally, the one standard deviation uncertainty in the 
strong coupling constant a s at Q 2 = M§ pOf contributes a 9% uncertainty at Q 2 = 30 
GeV 



2 



The variation of the parameters of the unfolding procedure within wide ranges had a 
negligible effect on the result. 

The Monte Carlo model dependence of the BGF and QCD-Compton acceptances 
was estimated on the hadron level by comparing the numbers obtained from the CDM, 
HERWIG and MEPS models and was found to be below 15% (see section 5). 

The matrix element cut-off in the Monte Carlo generator on the hadron level was 
varied by ±1 GeV in order to study how it affects the event composition. No significant 
change was observed. 

The radiative correction has been estimated for each x g / p bin by calculating the cross 
section with and without including QED radiation in the DJANGO 2.1 Monte Carlo pro- 
gram |^T | . This led to variations consistent with the statistical accuracy of the generated 



Monte Carlo samples. 

The full analysis was repeated using the JADE jet finding algorithm [22] in the labora- 



tory frame with a fixed mass cut of 10 GeV. This is an algorithm based on a fundamentally 
different principle of reconstructing jets compared to the cone algorithm which results in 
different numbers of selected events and different migrations. In spite of this the finally 
obtained results on the gluon density agree within the statistical accuracy. The reason 
for using the cone algorithm to extract the final results was the better x g / p resolution and 
the greater suppression of tails. 

As a further systematic check the analysis was repeated with a stricter cut in the 
invariant mass of the hard sub-system (s > 200GeV 2 ), again yielding variations well 
within the statistical accuracy, although the statistical errors were significantly increased. 



8. Discussion of the Results 

So far only indirect constraints on the gluon density based on a Q 2 evolution of the quark 
density have been obtained for low values of x. A comparison between such indirect 
extraction and direct measurement of the gluon density constitutes an important test of 
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Figure 5: The measured gluon density at an average Q 2 of 30 GeV 2 as a function of the 
fractional gluon momentum compared with indirect determinations by HI and ZEUS 



at Q 2 = 20 GeV 2 as well as with a determination from J/\l/ production by NMC [27] 
evolved to Q 2 = 30 GeV 2 (see text). 



perturbative QCD. Fig. 5 shows a comparison of this measurement with recent indirect 
determinations by the HI |23| and ZEUS J24j collaborations at a Q 2 of 20 GeV 2 which 
covers an even lower x g / p region. Our results are also consistent with a previous indirect 
measurement by NMC |25|] in a region x g / p > 10~ 2 . 



Direct measurements have been reported by UA2 [^] and NMC |?7| at x g / p values 
above 0.04. The NMC data based on inelastic J/ip production in DIS is included in Fig 5. 
The parametrisation of the NMC data provided by this collaboration was used to perform 
a LO DGLAP evolution from Q 2 = 1.5 GeV 2 to Q 2 = 30 GeV 2 . 

Further the results are compared with two different parametrisations of the gluon 
density in LO. The GRV model assumes the gluons and sea quarks to be valence-like at 
Q 2 = 0.23 GeV 2 and the growth of the gluon density with decreasing x g / p values is due to 
radiation of low x partons generated according to the Altarelli-Parisi evolution equations. 
The CTEQ3L parametrisations are based on input distribution functions at Q 2 = 4 GeV 2 
assuming the sea and gluon distributions to have the same power dependence of the 
fractional momentum (x Bg ). Both are consistent with our result. 

The contribution of diffractive events, in which the exchanged boson interacts with 
a colourless object in the proton and therefore does not span a colour string between 
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the hard sub-system and the proton, has been evaluated using data by selecting events 



with no activity in the forward detectors [pS j. After correcting for the efficiency of the 
diffractive selection and for non diffractive events passing the selection cuts, a (8 ± 2)% 
contribution is observed in the data sample consistent with a flat x g / p dependence (see 
also Table 1). 

Next to leading order (NLO) corrections using the cone algorithm for jet production 
in deep-inelastic scattering do not exist. In photoproduction these corrections have been 
found to be of the order 25% for an average < ]f T > (see Tab. 1) similar to this analysis 
|29fl . Although in photoproduction besides the contribution from BGF and QCDC the so 
called resolved processes have been taken into account, NLO corrections of similar size 
are expected 



9. Conclusions 



From a measurement of the cross section for (2+1) jet events, a direct LO determination 
of the gluon density in the proton has been performed in a previously inaccessible domain 
of the gluon fractional momentum x g / p (1.9 • 1CT 3 < x g / p < 0.18) at an average value of 
Q 2 = 30 GeV 2 . Our data are consistent with a steep rise in the gluon density as x g / p 
decreases. Recent indirect extractions of the gluon density based on the Q 2 evolution of 
the structure function Fi by the ZEUS and HI experiments are compatible with our data. 
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